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Abstract 

Starting from the divergence pattern of perturbative quantum chromody- 
namics, we propose a novel, non-power series replacing the standard expansion 
in powers of the renormalized coupling constant a. The coefficients of the new ex- 
pansion are calculable at each finite order from the Feynman diagrams, while the 
expansion functions, denoted as W n (a), are defined by analytic continuation in 
the Borel complex plane. The infrared ambiguity of perturbation theory is man- 
ifest in the prescription dependence of the W n (a). We prove that the functions 
W n {a) have branch point and essential singularities at the origin a = of the 
complex a-plane and their perturbative expansions in powers of a are divergent, 
while the expansion of the correlators in terms of the W n (a) set is convergent 
under quite loose conditions. 



1 Introduction 



The standard QCD perturbative expansion in powers of the renormalization group 
improved coupling is plagued by several deficiencies. The series is divergent and even 
Borel non-summable, having a zero convergence radius fl]]-[|| and coefficients exhibit- 
ing asymptotically factorial growth and nonalternating signs 0-||. The truncated, 
fixed order perturbative expansion is afflicted with the problem of renormalization 
scheme dependence and violates explicitly, due to the Landau singularities, the rigor- 
ous momentum plane analyticity imposed by general principles for confined theories 
0. Several resummations or reformulations of the perturbative expansion have been 
considered recently, trying to cure or reduce these deficiencies. Modified expansions 
which eliminate the unphysical singularities in the momentum plane were proposed in 
T0| , [|TTJ . Reformulations of the perturbation theory attempting to reduce the renor- 
malization scheme dependence of the truncated series were also proposed by several 
authors, using concepts like minimal sensitivity |T2|, effective charge [13], |14| or Pade 
approximants in the complex plane of the coupling constant |15[-||T7 . 

Other attempts to go beyond the conventional perturbative expansion are moti- 
vated by the large order behaviour of perturbation theory |f[8[| -[p2| . The perturbation 
theory is intrinsically ambiguous, due to the infrared regions of the Feynman diagrams. 
According to current views [^1], the perturbative ambiguity will be compensated in the 
complete theory by the nonperturbative contributions. But the estimate of this intrin- 
sic ambiguity is made difficult by the fact that the series is divergent, and the higher 
terms dramatically spoil the accuracy of the result above a certain order. If one suc- 
ceeded to replace the usual perturbative expansion by a convergent series (each new 
term improving the accuracy of the approximation instead of spoiling it), the intrinsic 
ambiguity of perturbation theory would be better defined. In the present work we 
address this problem and propose a new, non-power perturbative expansion in QCD, 
using the principle of analytic continuation in the Borel complex plane. 

Although the perturbative series in QCD is not Borel summable, since the condi- 
tions for Borel summability [23], [24] are not satisfied, the notion of Borel transform of 



the QCD correlators is nevertheless useful, as its singularities in the Borel plane con- 
tain much physical information. The infrared (IR) and ultraviolet (UV) renormalons 
and the instantons define a generic doubly-cut Borel plane, which provides an intuitive 
measure of the ambiguities of the perturbation theory. 

The conformal mapping of the Borel plane was suggested in [25] as a technique to 
reduce or eliminate the ambiguities (power corrections) due to the large momenta in the 
Feynman integrals. Such conformal mappings, which move the UV renormalons further 
from the origin f26 |, [ 27f exploit however only in part the known singularity structure 
in the Borel plane. An optimal mapping, which performs the analytic continuation 
in the whole doubly-cut Borel plane, was considered in our previous papers [p8[| , |29j 
(similar mappings were applied later in [30|, |3l]]). By means of this technique, we 
defined in [^] a non-power perturbative expansion in QCD in terms of a new set of 
functions that fully exploit the location of the singularities in the Borel plane. 

In the present paper we study the properties of these expansion functions, denoted 
below as W n (a), and the role of the modified expansion for a better understanding of 
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the ambiguities inherent in perturbative QCD. The results are impressive: in contrast 
to the perturbative powers a n , the new functions W n (a) provide a convergent expansion 
of the correlator in a large domain of the coupling constant complex plane, and share 
in addition certain important properties with the expanded correlator itself. 

The paper is organized as follows: in section 2 we recall the basic facts known about 
the singularities of the QCD correlation functions in the Borel plane and discuss the 
method of optimal conformal mapping. We use for illustration the Adler function in 
massless QCD, but the generalization to other cases is straightforward. In section 3 
we define the functions W n (a) and the modified perturbative expansion replacing the 
standard series in powers of the coupling constant. We investigate the analytic prop- 
erties of the expansion functions in the complex a-plane, their asymptotic expansions 
for small a, and the convergence conditions for the new expansion. In the concluding 
section 4, we point out the specific merits of the modified expansion for defining the 
intrinsic ambiguity of perturbation theory. 

2 Conformal mapping of the Borel plane 

2.1 Singularities in the Borel plane 

We consider the Adler function in massless QCD 

D = -4vr 2 s 




where H(s) is the amplitude of the electromagnetic current-current correlation function 

UT(p) = i J d 4 xe" ipx, (0| T{f{x)f{0)) |0) 

= (g^p 2 -p^Ms), s = p 2 . (2) 

In perturbative QCD, the function D can be formally expanded in powers of the renor- 
malized coupling constant a s {ii 2 ) 

(3) 

where the coefficients D n depend on the external momentum squared s, and are renor- 
malization scheme and scale dependent. Certain classes of Feynman diagrams suggest 
that they are factorially increasing with n, therefore the series is divergent, and has 
to be given a precise meaning. 

Among various summation techniques of power series [23], [54], the Borel method 




received much interest in recent years, although the mathematical conditions required 
for its use are not satisfied in QCD Jl]], [[J. As mentioned in the Introduction, even 
if the Borel summation method is not applicable it is useful to examine and exploit 
the singularity structure of the Borel transform, as it contains important physical 
information. In the present paper we use the Borel transform to define an improved 
perturbative expansion, using in addition the principle of analytic continuation. 
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The Borel transform of the Adler function is defined by the power series 



B(u,s) = J2b n u n } (4) 

n=0 

with b n related to the original perturbative coefficients appearing in (|3]) by 

where (3q is the first coefficient of the (3 function. Then the Adler function given by the 
series @ can be formally expressed as the Laplace transform 

oo 

D(s, a) = - [ e~ u/a B(u, s) du , (6) 
a J 
o 

where a = Poa s (fj, 2 )/n. 

The Borel transform B(u,s) is scheme and scale dependent. As discussed in [2T|, 
in the large Nf limit the momentum dependence can be factorized as 

Cc \ ~ u 
B o(u) , (7) 

where C is a scheme-dependent constant. The power factor in (|7|) combines with the 
exponential in ([]) in a renormalization scheme and scale independent quantity, the 
remaining factor Bq(u) also being scheme and scale independent. The validity of a 
similar factorization beyond the large Nf limit is an open question [21]. However, it 



is generally assumed that the position of the singularities of the function B(u,s) in 
the Borel u plane is independent of renormalization scheme and scale, and also of the 



external momenta (as discussed in |2T[], this is expected to be the case at least in the 
so-called "regular schemes"). 

The method developed in the present paper is based on the analytic continuation in 
the Borel plane, requiring therefore the knowledge of the singularity structure in this 
plane. We shall assume that the location of the singularities of the Borel transform is 
known and is independent of the renormalization scheme. For the Adler function, the 
coefficients D n are assumed to have the specific large order increase 



'■ -xo..^ i^y . (8) 

k 



where the k in the sum are integers, with and depending in general on s and /z 2 . 
This behaviour leads to branch point singularities for the function B in the w-plane, 
along the negative axis (the ultraviolet renormalons) and the positive axis (the infrared 
renormalons) . Specifically, the branch cuts are situated along the rays u < — 1 and 
u > 2 (see Fig. []]), and the nature of the first branch points was established in [25] 



and in |32]| . The additional singularities due to the instanton-anti-instanton pairs [33[] 
are situated at larger positive values of u, and will not influence the method discussed 
in this paper. 
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Figure 1: The Borel plane for the Adler function. 
2.2 Optimal conformal mapping 

The singularities of the Borel transform produce the perturbative ambiguities of the 
QCD correlators. For the Adler function, the Taylor expansion (f|) of the Borel trans- 
form converges on the disk \u\ < 1 that reaches the nearest singularity, which is the 
first UV renormalon. The corresponding power correction ambiguity is however related 
to large momenta in the Feynman diagrams, and can be minimized or eliminated in 
QCD within perturbation theory. As suggested in Ref. |25], this can be achieved by an 



analytic continuation outside the convergence circle, by means of a conformal mapping. 



This method was investigated in ^7|, using the mapping 



y/l+U-1 , 

Vl + u + 1 

and the corresponding expansion 

oo 

B(u,s) = Y,d n v n (u), (10) 

n=0 

instead of the power series (||). The mapping @ transforms the UV cut u < — 1 in 
the w-plane onto the unit circle \v\ = 1 in the v plane, while the IR cut u > 2 remains 
inside the disk \v\ < 1 (see Fig. |^). As a consequence, the series fllCf) converges inside 
the small circle, which touches the image of the lowest IR renormalon in the f-plane. 
Thus, although the UV ambiguities (power corrections) are in principle eliminated by 
this method, the series flT0| ) remains divergent along the positive semiaxis, the IR cut 
being outside the convergence circle. 

The analytic continuation of the perturbative Borel transform into the whole holo- 
morphy domain D can be performed by an optimal conformal mapping, w(u), which 



maps the whole D onto the unit disk [34|. For the Adler function this mapping is [28 



y/T+u- Jl-u/2 
w(u) = 1=7= • 11 
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Figure 2: The f-complex plane. The UV cut is mapped upon the unit circle, while the 
IR cut is situated inside it. The convergence domain of the series fliTf ) is limited by 
u = 2, the lowest branch point of the IR cuts family. 

The two branch points of w(u), u — — 1 and u = 2, coincide with the lowest branch 
points of the Borel transform B(u,s), while the corresponding two cuts, u < — 1 and 
u > 2, cover the other branch points of B. As shown in Fig. || these two rays are 
mapped by (|TT| ) onto the boundary circle of the unit disk \w\ = 1 in the w— plane. The 
inverse of the mapping (0) is 



3w 2 -2w + 3 3(w-()(w-(*) ' 1 ' 

where ( = (\/2 + %)j (y/2 — i) and its complex conjugate (* are the images of u = oo 
on the unit circle. 

As proposed in [ ^Sfl , we expand B in powers of the variable w 

oo 

B(u,s) = Y / c n w n (u), (13) 

ra=0 

where the coefficients c n can be obtained from the coefficients bk, k < n, using Eqs. (f|) 
and (|TTD . By expanding -B(m, s) according to (p!3|) one makes full use of its holomorphy 
domain, because the known part of it (i.e. the first Riemannian sheet) is mapped onto 
the convergence diskQ The series (|13]) converges inside the whole disk \w\ < 1, i.e., 
in the whole cut u plane, up to the cuts. A very important property, proved in [34], 
is that the expansion ([i~3|) has the fastest asymptotic (large-order) convergence rate, 
compared to any other expansion in powers of a variable that maps only a smaller 

lr The conformal mapping used here is confined to the first sheet. If the branch point structure and 
analyticity on other sheets were known, a multi-sheet mapping could be used which simultaneously 
uniformizes several branch points |3q|. 
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Figure 3: The w;-complex plane. Both the UV cut and the IR cut are mapped upon 
the boundary unit circle. Inside it, there is no singularity. The convergence domain is 
the whole (cut) Borel plane. 



part of the holomorphy domain onto the unit disk. We recall that the large-order 
convergence rate of a power series is equal to that of the geometrical series with the 
quotient r/R, r being the distance of the point from the origin, and R the convergence 
radius. The proof given in [34j| consists in comparing the magnitudes of the ratio r/R 



for a certain point in different complex planes, corresponding to different conformal 
mappings. When the whole analyticity domain T> of the function is mapped on a disk, 
the value of r / R is minimal [[!4[ . 

3 Modified QCD perturbative expansion 
3.1 Non-power expansion functions in QCD 

The expansion flT3] ) of the Borel transform suggests to expand the Adler function in 



the series [2S, 



D(s,a) = Y,CnW n (a) } (14) 

n=0 

where the functions W n (a) are defined as Laplace transforms of w n (u) : 

oo 

WJa) = - I e~ u/a w n (u) du, n = 0, 1, 2, . . . . (15) 
a J 



At each finite truncation order N, the expansion (II) is obtained by inserting the 



series (|13|) into the Laplace integral (§) and exchanging the order of summation and 
integration. This procedure is trivially allowed at any finite integer N > 0. For N — > 
oo, however, the new expansion ( fUl) represents a nontrivial step out of perturbation 
theory, replacing the perturbative powers a n by the functions W n (a). As proved in |2Df , 
the expansion ( |14|) converges under certain rather general conditions (we discuss this 
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point below in subsection |3.4|) . The convergence of the series (|I4D is our key argument 
in favor of the stepping out of the standard perturbation theory. 

Our procedure is an obvious generalization of the conformal mapping method pro- 
posed in [3H] for Borel summable functions. Formally, the expansion (H) is obtained 



from the standard perturbative expansion (0) by replacing the coefficients b n , appearing 
in the Taylor series (§), by the coefficients c n of the improved expansion (|13D , and the 
perturbative functions n\a n (which multiply the coefficients h n ) by the new functions 
W n (a) defined by the integral (|l"5l) . Actually, this integral is not well-defined since the 
branch-point u = 2 is encountered along the integration range. This is a manifestation 
of the intrinsic ambiguity of the perturbation theory produced by the infrared renor- 
malons. In defining the functions W n , we shall use the same prescription as the one 
adopted in ([]) for the correlator D itself. We shall consider in particular the functions 

W^\a) = - I e~ u/a w n {u) du , n = 0,1,2, ... , (16) 
a J 

c± 

where C + (C_) are lines parallel to the real positive axis, slightly above (below) it, and 
the principal value (PV) prescription 

Wr\a) = \[W^\a) + Wt\a)}. (17) 

In what follows we shall examine the expansion functions W n (a), showing that in many 
respects they resemble the expanded function D(s,a) itself. 



3.2 Analytic properties of W n (a) in the complex a-plane 

The problem of the analytic properties of the QCD correlators in the coupling constant 
plane is very complex, 't Hooft and Khuri || showed that renormalization group 
invariance and the multiparticle branch-points on the timelike axis of the momentum 
plane s = p 2 imply a complicated accumulation of singularities near the point a = 
0. Since the proof uses a nonperturbative argument (multiparticle states generated 
by confinement in massless QCD), it is difficult to see this feature in perturbation 
theory, even in partial resummations that take into account infinite classes of Feynman 
diagrams. 

In the present subsection we shall discuss the analytic properties of the expan- 
sion functions W n (a) defined by the integrals flip]). Since the integrand is bounded, 
|w n (-u)| < 1, the integrals converge for Rea > 0, i.e. —n/2 + 5 <ip < n/2 — 5, for 5 > 
arbitrarily small, where ip is the phase of a (a = \a\e 1 ^). Let us consider the functions 
W^ + \a) and take first a complex in the first quadrant, i.e. < ip < tt/2. Then we 
can rotate the integration contour C + in ( jl6|) towards the upper quadrant by an angle 
if), since the integrand is analytic. More exactly, we apply the Cauchy theorem for a 
closed contour going along C + up to \u\ = R, continued with the sector of the circle 
u = Re %< ^, with < <p < ip, and the segment u = re*^, with < r < R. This gives 

w n {Re irt> )d<P = 



1 - u / a w n (u)du+ — 



a 



a 



e l<t> exp 



-Re^-y) 
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R 



R/\a\ 



— / P -^/t a U/, n 



w re 



dr 



,— t „,,n 



w n (ta) dt. 



The second integral on the first line can be estimated as follows 



R_ 

\a\ 



cxp 



_Rcos(0 — iji) 



dcj) < ip 



R 


Rcosip 


- — - exp 


a 


a 



(19) 



because |w(«)| < 1 and < <fi < ip. Since ip < tt/2, the expression ( |T9"D vanishes for 
R — > oo. Taking this limit in ([18]) one obtains the equality 



"*w ra (to)dt, 0<V<vr/2. 



(20) 



Let us take now a in the fourth quadrant, i.e. —tt/2 < ip < 0. A representation of 
the form ( |20| ) can be obtained by making a rotation of the integration contour into the 
fourth quadrant up to the negative angle ip. But this rotation is not allowed for the 
contour C+, since one encounters the branch cut of the integrand w n (u). By Cauchy 
theorem it is however easy to relate W^ + ' to the function W^~\ defined by an integral 
along the contour C_. Taking into account the discontinuity of the function w(u) 
defined in (0), we have 



2i 



Wi +) {a) = WtHa) + — / e" u/a Im[w»] du 



(21) 



where 



\m\w n {u) 



4 _|_ n, \ n 1 2 J 

4 e (-i)^r 1 

p=0 



3m 



(l + tz)(u/2-l) 



l + u/A 



2p+l 



(22) 



is the imaginary part for u > 2 when u approaches the cut from above. It is convenient 
to make the change of variable u = 2 + y in the last integral of (pl|), writing it as 



2i- 



2 OO 

e 'a 



- y,a f n {y) dy. 



(23) 



where we denoted 

f n (y) 



• i + y/6 V 



E 

p=0 



•1) P C^ 



(3 + y)y/2 
3/2 + y/4 



2p+l 



(24) 



From (|22|) it follows that the f n (y) are algebraic functions, with poles and branch points 
for y < and cuts which can be taken along y < 0, so the integral (|23|) is defined. 
Moreover, we can rotate the integration axis to the lower quadrant, up to the angle 



S 



ip < 0. This rotation is allowed also for the function leading to a representation 

identical with the r.h.s. of (120). Collecting all the terms we obtain the expression 



CXJ OO 

Wi + \a) = J e- t w n (ta)dt +2ie-^y e~ l f n (ta) dt , -vr/2<^<0. (25) 
o o 

We can use similar arguments for the functions W^: the contour C_ can be rotated 
towards the fourth quadrant for negative ip, giving 

oo 

Wt\a) = J e~* w n (ta) dt , -tt/2 < tj; < , (26) 
o 

while for ip > one must first cross the real axis, picking up the contribution of the 
discontinuity of the integrand. This gives 



WH(a) = / e- t w n (ta)dt -2ie~ / e~' f n (ta) dt , < ip < tt/2 . (27) 
o o 

For the PV prescription defined in (|T7|), we obtain from Eqs. (p0[ ) and (^7|): 

oo oo 

W^ py )(a) = / e~*u; n (ta)dt -ie~» / e~*/ n (ta)dt, < t/j < vr/2 , (28) 
o o 

and from Eqs. fl2"5| ) and (^Bj): 

oo oo 

W!f v \a)= [ e- t w n (ta)dt +ie~l f f n (ta) dt , -tt/2 < ip < . (29) 



The expressions fl2"0|) and (|25| ) - ( |29|) define holomorphic functions in the right half 
plane Rea > 0, outside the real positive axis. When the coupling approaches this axis 
from above, i.e., for a + ie, with a > and e > 0, the expression (|28| ) becomes 



DO 



W^ py) (a + ^)= /" e -*ReK(ta)]dt, (30) 



since, by the definition (f24|), the imaginary part of the first integral is exactly cancelled 
by the last term in (p8|). Similarly, for a — ie the expression (|29D gives 



oo 



W { n PV \a-ie)= J e-*ReK(ta)]dt, (31) 
o 

where we used the reality condition w(u) = w*(u*) obvious from (|H]). Thus, the 
functions Wj^^ (a) have no discontinuity along the positive axis, where they take real 
values. 
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The expressions ( g8j ) and ( [29]) can be analytically continued from the right half plane 
to the left one, Rea < 0, up to the negative real axis, since the integrals converge and 
the integrands have no singularities along the integration contours. When a approaches 
the negative real axis, the imaginary part of the first term in Eq. Q28D (or (53)) is given 
by the cut of w(u) along the negative axis, and it is no longer cancelled by the imaginary 
part of the last term, given by the expression (|24]) analytically continued to negative 
values of the argument. So, the functions W^ PV \a) have branch cuts for a < 0. The 
real analyticity of the function w(u) implies however the equality 

Wr\a + ie) = (Wi PV \a-te)y. (32) 

Therefore, the W^ PV \a) are real analytic functions in the whole complex a-plane, 
except for a cut along the real negative axis and an essential singularity at a = 0, seen 
explicitly in the representations (|28|) and (p^). 

For other prescriptions the reality conditions written above are not satisfied. Thus, 
approaching the positive real axis from above and from below in (p0|) and (p5|), respec- 
tively, we obtain 

oo 

Wi +) (a ± ie) = J e - * {Re[w n (ta)] + i lm[w n (ta)]} dt , (33) 
o 

with Im[w n (£)] defined in fl22"|) . The functions W( + '(a) have no discontinuity along the 
positive axis, but they are complex for positive values of the coupling. They can also 
be extended into the whole complex plane, except for the essential singularity at a = 
and a branch cut along the negative axis. 

At each finite order of truncation in (|I~4D, the expanded function 



N 

I 

n=0 



DW(s,a) = J2c n W n (a) (34) 



will have the same analyticity properties as the W n (a). Thus, in the PV prescription, 
the D^ N \s, a) have a branch cut along the negative real axis and an essential singularity 
at the origin a = 0. It is instructive to see what are the consequences for the momentum 
plane analyticity, taking the renormalization point fi 2 = — s > and the one loop 
expression of the running coupling a(s) = l/ln(— s/A 2 ). With the usual definition of 
the logarithm on the first sheet (ln(— s) > for s < 0), the positive real axis a > 
corresponds to the part s < —A 2 of the space-like axis in the s-plane. We obtained 
no branch cuts in this region, in agreement with the requirements of unitarity and 
causality for QCD On the other hand, the negative real axis a < corresponds 
to the Landau region —A 2 < s < 0. The W n (a) have a branch cut along this region, 
obtained explicitly by analytic continuation using the representations ( p8|) and ([291). 
So, in finite orders, the new expansion (|I~4T), calculated outside the Landau region, can 
be analytically continued inside the Landau region, where it exhibits a branch cut. 



Similar properties were obtained in the large (3q limit in p? 



10 



Figure 4: The functions W m (a) defined with the PV prescription fllTD , m = 1, . . .6 
for a from to 0.8. Unlabelled are the perturbative expansion functions n\a n for 
n — 1, ... ,6. 



3.3 Asymptotic expansion of W n (a) for small a 

In this subsection we investigate the perturbative expansion of the functions W n (a) in 
powers of a. As shown in the previous subsection, W n (a) have singularities at a = 0, 
so their Taylor expansions around the origin will be divergent series. We take first a 
real and positive. The asymptotic expansion is obtained by applying Watson's lemma 
38fl . Specifically, we consider the Taylor expansion 



w n (u) = £ tfV , (35) 

k=n 

which is convergent for \u\ < 1. The sum begins with k = n since, as follows from flllD , 
the derivatives (w n )( k \0) vanish for k < n (in particular = (3/8)™). Choosing a 
positive number X < 1, we can express w n (u) for < u < X as 

JV 

^» = Eei n V+^, (se) 

k=n 

with a bounded remainder 

?( n )| ^ T\/T 



R^\<M n u N+ \ (37) 



We write now W n (a) as 



X , oo 



W n (a) = - [ e~ u/a w n {u) du+ - [ e" u/a w n (u) du , (3* 

X 



and insert in the first term the expansion (|36|) . This gives 



X , , X oo 



W n a = - 
a 



- / e" u/a V^Vd«+^ / e-"/^ nM 7V+1 d M +- /e- M / a W "( M )d M , (39) 
a J , a J a J 

n k=n n v 
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where \9 n \ < 1. Since |w n (u)| < 1, the last term in (|39|) is bounded as 



l/a f e~ u/a w n (u)du 



< e 



-X/a 



(40) 



On the other hand, for fixed k > we have 

x 



- [ e~ u/a u k du = k\a k --( e~ u/a u k du . 
a J a J 

o x 



(41) 



In the last term we make the change of variable u = X(l + y) and use the inequality 
(1 + y) k < e ky for k > [38 , which implies that the following estimates 



x 



e -u/a u k du 
X k+i 



X k+i 



a 



X/a / y(X/a-k) 



e -X/a J e -Xv/*(l +y )k dy < 


X k+1 



dy 



X/a 



a(X/a — k) 

u 

are valid for fixed k and small a. Therefore Eq. (|41| ) can be written as 



(42) 



x 

- J e- u/a u k du =k\a k + O (e~^) . 



(43) 



By using this estimate in ( |39"|) and taking into account the inequality (f40l) , we obtain 



N 

W n (a) = J] i { k n) k\a k +M n (N + l)\ a N+1 + O (e~£) , 
where M n is independent of N. One has therefore the asymptotic series 



(44) 



(")iuufc 



a -> 0, . 



(45) 



k=n 



It is easy to extend the above arguments for complex values of a in the right half-plane, 
\i>\ < n/2 -5 with 5 > 0. 

The expansion (53) is independent of the prescription required in the definition 
of W n . (It is formally obtained by inserting the expansion (|35|) in Eq. (pl|) and 
integrating term by term from to oo.) We notice that the first term of each W n (a) 
is proportional to n\a n with a positive coefficient, thereby retaining a fundamental 
property of perturbation theory But the series (^3J) are divergent: indeed, since the 
expansions (^) have their convergence radii equal to 1, then for any R > 1 there are 
infinitely many k such that > R~ k ||38|| . Actually, the divergence of the series ( ^5|) 
is not surprising, in view of the analyticity properties derived in the previous section. 
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Figure 5: The same graph as in Fig. |J with a in the interval (0, 6). 

For illustration we give below the expansions of the first W n (a) 

Wi(a) ~ 0.375a - 0. 187a 2 + 0.457a 3 - 1.08a 4 + 4.32a 5 + ..., 
W 2 (a) ~ 0.281a 2 - 0.422a 3 + 1.58a 4 - 5.80a 5 + 29.78a 6 + ... , 
W 3 {a) ~ 0.316a 3 - 0.949a 4 + 5.04a 5 - 25.95a 6 + 167.99a 7 + ... . (46) 

The higher powers of a become quickly important in (pEB]), the expansion coefficients 



eventually adopting factorial growth. For instance, the 5th-order coefficients in (f46|) 
all equal 5 approximately, while the lOth-order ones are between 5 x 10 4 and 9 x 10 4 , 
with alternating signs. The functions W n (a) have divergent perturbative expansions, 
resembling the expanded QCD correlation function D. 

Although the series ( f4~6"D are divergent, the functions W n (a) are well-defined (once 
a prescription has been adopted), and bounded in the right half plane Rea > 0: 

|W n (a)| < ^ re- uRea ^ 2 \w n (u)\du < P~ , (47) 
\a\ Jo Kea 

since |w n (M)| < 1. For a real and positive the right hand side of ([47]) is equal to unity. 
We indicate in Figs. f| and |5] the shape of the first functions W n , calculated with the 
PV prescription, for real values of a. 



3.4 Large-order behaviour of W n (a) and convergence condi- 
tions 

The convergence of the series ([TJ|) is not a priori obvious. Indeed, the expansion (O) of 
the Borel transform in powers of w converges for points \w\ < 1 arbitrarily close to the 
integration axis, but not necessarily on the boundary. Intuitively, one may expect that 
the convergence depends on the strength of the singularities of the Borel transform. 

We investigated the convergence problem in ||29|| , by estimating the behaviour of 
W n (a) for large n with the method of steepest descent |36|], [[3811 . The function W^ + \a) 
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defined in 



can be written as 

W&(a) 



F »V> du. 



c± 



where 



i i u 

F n {u) = ~ 

a 



n Ytiwyu) . 

The saddle points of the integrand are given by the equation 

w'{u) 1 



w[u) 



an 



(4£ 



(49) 



(50) 



which has four complex solutions. The technique applied in consists in rotating 
the integration axis in the complex u-plane, without crossing singularities, up to the 
nearest saddle points (located in the first (fourth) quadrant for C + ( C_ ), respectively). 
Omitting the details given in PD|] , we quote the asymptotic behaviour of W^ + '(a) for 
n — * oo 



i/-n p -2 3 / 4 (l+i)(n/a) 1 /2 



ni( n e 



n 



oo . 



(51) 



with ( defined below (|T|). Similarly, the large n behaviour of \a) is 



*\n -2 3 / 4 (l-i)(n/a) 1 / 2 



These expressions are valid for complex a, a = 

|V>] < 7T/6. 

The convergence of the expansion ( ]T4] ) depends on the ratio 

c„^ ±} (a) 



n — > oo . (52) 
e 1 ^, where ifi satisfies the inequality 



(±) 



(53) 



(54) 



As discussed in P5| , if the coefficients c n satisfy the condition 

|c n | < Ce enl/2 (55) 

for any e > 0, then the expansion fll4|) converges for a complex in the domain 

Re[(l ± i)cr 1/2 ] > 0, (56) 

which is equivalent to < tt/2 — 5. Since the condition ( j53|) is more restrictive, it 
follows that, if the condition ( p3| ) is satisfied, the series (|T4T ) converges in the sector 
defined by (|53|). 

The coefficients c n are obtained by inserting into the Taylor series (|]) the expan- 
sions in powers of w of the function u(w) defined in ([12]). A precise estimate of the 
behaviour of the c n starting from a general form of the standard perturbative coeffi- 
cients D n is difficult to obtain. In the special case of a Borel transform with branch 
point singularities, considered in | 29| , one can derive the generic behaviour 

\c n \ < C'n^ = C'e^ lnn , £>0, (57) 
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which satisfies the convergence condition flog) . Whether this bound is valid or not in 
general in QCD is an open problem, but Q57|) nevertheless represents a rather conser- 
vative assumption. 

In subsection we established analyticity properties for the functions W n (a) in 
the complex a plane. The expanded function D(s,a) will be holomorphic in the do- 
main where the series (|14]) converges uniformly. Using the results of subsection |3.2| , 
we conclude that D(s,a) is expected to be holomorphic inside the sector (|53| ) of the 
complex a-plane, with an essential singularity at the origin. This region is larger than 
the horn shaped domain found in PJ and [fj] . The difference can be easily understood: 
indeed, the multiparticle singularities of the correlators in the momentum plane for 
s — > oo, essential for the horn shaped boundaries preventing the analytic continuation 
near a = 0, ||, cannot appear in perturbative massless QCD, which is the frame 
adopted here. 



4 Concluding remarks 

In the present paper we used the analytic continuation in the Borel plane to define a 
modified perturbative expansion in QCD, the n-th order coefficient c n being calculable 
from the standard perturbative coefficients D^, k < n. The new expansion functions 
W n (a) replacing the powers of the coupling constant are defined in ( |I5|) as the Laplace 
transforms of the n-th power of the function w(u), which performs the conformal 
mapping of the whole doubly-cut Borel plane onto a unit disk. 

The definition of the integral (|15|) is ambiguous due to the branch cut u > 2 along 
the integration path. We propose to define W n (a) by choosing the same prescription 
as that used to define the expanded function D(s) itself, (||), so as to make the W n (a) 
resemble D(s) as much as possible. The form of the W n (a) is not accidental or chosen 
ad hoc: it is imposed by the analyticity properties of the Borel transform, which in 
turn are determined by the information contained in the perturbative coefficients of 
the QCD correlation functions at all orders. 

We proved that 

/i/ the functions W n (a) are analytic in the complex a-plane, with essential and 
branch point singularities at a = (subsection 3.2), 

/ii/ their expansions in powers of the coupling constant a are asymptotic divergent 
series, the lowest order term of W n (a) being proportional to a n (subsection 3.3), and 

/iii/ the expansion (|T3[) of the correlation function in terms of the W n (a) set con- 
verges under plausible conditions, such that are expected to be satisfied in perturbative 
QCD (subsection 3.4, see also our previous paper |2"9[] ). By this we avoid the fatal di- 
vergence of the standard perturbative expansion: for the new expansion in terms of 
the W n (a), the addition of higher-order terms does not damage the result, as is the 
case with the series in powers of a. The convergence property is the main merit of the 
new expansion. 

The prescription required for calculating the Laplace integral in the definition ( jI5|) 
reflects the intrinsic ambiguity of perturbation theory, which originates from the in- 
frared regions of the Feynman diagrams and is manifest in the presence of the IR 
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renormalon cut. Once a definition of the perturbative ambiguity is adopted (as the dif- 
ference between two integration prescriptions, for instance), the series ( |HD provides a 
systematic calculation of this ambiguity at higher orders, since the expansion in terms 
of the W n (a) is convergent. This is to be contrasted with the standard expansion, 
where the ambiguity is defined in a less precise way by truncating the series at some 
finite order, beyond which the terms start to increase. Unlike in this procedure, the 
effects of the UV and IR parts of the Feynman diagrams have been disentangled in our 
approach. 

Thus, using the analytic continuation in the Borel plane, we have been able to sep- 
arate two problems that are usually interconnected: the divergence of the series (which 
can be solved within perturbation theory), and the problem of the intrinsic infrared 
ambiguity of perturbation theory. This ambiguity, expressed in the prescription de- 
pendence of the correlator and the expansion functions W n (a), can be removed only 
when nonperturbative effects are included. 
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